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Abstract We briefly discuss the collinear factorization formula for the associated produc- 
tion of one particle and a Drell-Yan pair in hadronic collisions. We outline possible applica- 
tions of the results to three different research areas. 
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1 Introduction 

The description of particle production in hadronic collisions is interesting and challenging 
in many aspects. Perturbation theory can be applied whenever a sufficiently hard scale char- 
acterizes the scattering process. The comparison of early LHC charged particle spectra with 
next-to-leading order perturbative QCD predictions 1 1] shows that the theory offers a rather 
good description of data down to hadronic transverse momentum of the order of a few GeV. 
For hadron production at even lower transverse momentum, the theoretical description in 
terms of perturbative QCD breaks down since both the coupling and partonic matrix ele- 
ments diverge as the transverse momentum of final state parton vanishes. In this paper we 
will study the semi-inclusive version of the Drell-Yan process, H\ + H2 — > H + 7* + X, in 
which one particle is measured in the final state together with the Drell-Yan pair. In such 
a process the high invariant mass of the lepton pair, Q 2 , constitutes the perturbative trigger 
which guarantees the applicability of perturbative QCD. The detected hadron H could then 
be used, without any phase space restriction, as a local probe to investigate particle produc- 
tion mechanisms. The evaluation (2]|3] of ff{a s ) corrections shows that there exists a class of 
collinear singularities escaping the usual renormalization procedure which amounts to reab- 
sorbing collinear divergences into a redefinition of bare parton and fragmentation functions. 
Such singularities are likely to appear in every fixed order calculation in the same kinemat- 
ical limits spoiling the convergence of the perturbative series. In Refs 15][3) a generalized 
procedure for the factorization of such additional collinear singularities is proposed. The lat- 
ter is the same as the one proposed in Deep Inelastic Scattering [4| where the same collinear 
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singularities pattern is also found, confirming the universality of collinear radiation between 
different hard processes. Such a generalized factorization will make use the concept of frac- 
ture functions and the corresponding renormalization group equations |5 1. In a pure parton 
model approach, these non-perturbative distributions effectively describe the hadronization 
of the spectators system in hadron-induced reactions. The resulting transverse-momentum 
integrated cross-section is finite and valid for all transverse momentum of the detected 
hadron, without any restriction imposed by the singular behaviour of matrix elements. We 
further note that this process is the single-particle counterpart of electroweak-boson plus 
jets associated production |6|, presently calculated at nex-to-leading order accuracy with up 
to three jets in the final state | 7 1. One virtue of jet requirement is that it indeed avoids the 
introduction of fragmentation functions to model the final state, which are instead one of 
the basic ingredients entering our formalism. At variance with our case, however, jet recon- 
struction at very low transverse momentum starts to be challenging |8i] and it makes difficult 
the study of this interesting portion of the produced particle spectrum. 



2 Collinear facrotization formula 

The associated production of a particle and a Drell-Yan pair in term of partonic degrees of 
freedom starts at &{(X S ). One of the contributing diagrams is depicted in Fig. 10. Assuming 




Fig. 1 Example of diagram contributing to hadron production in the central fragmentation region to order 
ff{a s ) in eq.fl}. 



that the hadronic cross-sections admit a factorization in term of long distance non-perturbative 
distributions and short distance perturbative calculable matrix elements for the partonic pro- 
cess i(pi) + j{pz) l{k) + 7*> predictions based on perturbative QCD are obtained con- 
voluting the relevant partonic sub-process cross-sections, dd 1 ^ 1 ' , with parton distribution 
functions, fi and fj, and fragmentation functions, D H l l . The hadronic cross-sections, at cen- 
tre of mass energy squared S, can be symbolically written as [ 2 , 3 1 



where the convolution are over the momentum fractions of the incoming and outgoing par- 
tons. The variable z is the energy of the observed hadron H scaled down by the beam energy, 
VS/2, in the hadronic centre of mass system and p is its partonic analogue. The invariant 
mass of the virtual photon is indicated by Q 2 . The partonic indeces i, j and / in the sum 
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Fig. 2 Parton model formula, eq.{2)> for the associated production of a particle and a Drell-Yan pair. 

run on the available partonic sub-process. The superscripts label the incoming hadrons and 
the presence of crossed term is understood. Within this production mechanism, the observed 
hadron H is generated by the fragmentation of the final state parton /, and for this reason we 
address it as central. The amplitudes squared [9|, however, are singular when the transverse 
momentum of the final state parton vanishes. In such configurations, the parent parton I of 
the observed hadron H is collinear the incoming parton i or /. Perturbation theory looses 
its predictivity as these phase space region are approached. The same pattern of collinear 
singularities are found also in an analogue calculation in Deep Inelastic Scattering |'4J. In 
both processes such singularities can not be treated with the usual renormalization proce- 
dure which amounts to reabsorb collinear divergences into a redefinition of bare parton and 
fragmentation functions. A widely used procedure to deal with such singularities (for exam- 
ple in modern Monte Carlo generators) is to introduce an arbitrary cut-off on the produced 
parton transverse momentum. It would be, however, highly desiderable to develop a tech- 
nique to resum such logarithmic contributions to all orders in perturbation theory since such 
divergences will occur in every higher order calculation spoiling the convergence of the 
perturbative expansion in this kinematical limit. Fracture functions together with their own 
renormalization group equations (3J can be shown to provide the correct tool to perform 
such resummation. Bare fracture functions, Mf' Hl (x,z), parametrize hadron production at 
vanishing transverse momentum. They express the conditional probability to find a parton i 
entering the hard scattering while an hadron H is produced with fractional momentum z in 
the target fragmentation region of the incoming hadron. The use of fracture functions opens 
the possibility to have particle production already to ff((X® ), since the hadron H can be 
non-pertubatively produced by a fracture function M itself. The lowest order parton model 
formula can be symbolically written as 



da H,T,(0) r d rdXZr.Jl], ,J2], s , „[2], w[l], .ndd'j^ 

where the superscripts indicate from which incoming hadron, Hi or Hi, the outgoing hadron 
H is produced through a fracture functions. This production mechanism is sketched in 
Fig. ©. So far we have only considered &{a s ) corrections in the central fragmentation 
region, eq. (fl}, to the parton model formula, eq. (O . In order to complete the calculation 
to ff(a s ) we should also consider higher order corrections in process initiated by a fracture 
functions. In this case, in fact, the hadron H is already produced by these distributions and 
therefore final state parton in real emission diagrams should be integrated over and results 
added to virtual corrections. One of the contributing diagrams is depicted in Fig. The 



Fig. 3 Example of diagram contributing to ff(a s ) corrections in the target fragmentation region, eq.(5J. 



general structure of these terms is 

da H,T,(\) r dx r dx { ] { ] ^ AM , M^' 1 ' m 

-agar •*£/ — 7 — N W)/J 'fcO+«i W)/J ® 

We refer to this corrections term as to the target fragmentation contribution. The calcula- 
tion is, a part from different kinematics, completely analogue to the inclusive Drell-Yan 
case. The factorization procedure is accounted for by substituting in eq. ® the bare frac- 
ture and distributions functions by their renormalized version 03). Renormalized partem 
distributions and fracture functions homogeneous terms do cancel all singularities present 
in eq. {3}. The additional singularities in eq. (|T]l are cancelled by the combination of parton 
distributions and fracture functions inhomogeneous renormalization terms. Adding all the 
various contributions, the resulting /?, -integrated cross-sections, up to order 0(oc s ), is then 
infrared finite 02) and can be simbolically written as 

|pL « aoE [Mfl ®ff + (1 2)] (1 + gc'V) + g a ZfP ®ff SD*/^, 

(4) 

where do = Ana 2 m /9SQ 2 . We refer to the previous equation as to the collinear factorization 
formula for the process under study. The next-to-leading coefficients C' J and K l J have been 
calculated [3], making the whole calculation ready for numerical implementation. In the 
following section we outline some applications of the proposed formalism. 



3 Section Applications 

The first possible benchmark process for the calculation could be strangeness production 
associated with a Drell-Yan pair, p + p — >V + y*+X, where V generically indicates a A 
or A hyperon. The first process, according to so called leading particle effect, should be 
sensitive to the spectator system fragmentation into A hyperon at very low transverse mo- 
mentum and therefore to be modelled with the help of fracture functions in da H > T . In the 
anti-hyperon production case no such effect should be present so that A are expected to 
be mainly produced by the fragmentation of final state parton as described by da H C . The 
measurament of such process for different Q 2 could allow to test the fracture functions scale 
dependence embodied in their peculiar evolutions equations and the validity of the pro- 
posed factorization formula in all its components. With this respect, a comparison of strange 
particle production in hadronic collisions and DIS would be extremely interesting since lon- 
gitudinal momentum spectra of A and A hyperon has been measured quite accurately in a 



5 



number of charged and neutral current DIS experiments, both in the current as well in the 
target fragmentation region. 

The formalism may found application in the study of single diffractive hard process, 
p + p— > p + Y* + X, where the outgoing proton has almost the incoming proton energy and 
extremely low transverse momentum with respect to the collision axis. This process has been 
intensively analyzed in the DIS at HERA, revealing its leading twist nature. From scaling 
violations of the diffractive structure functions 1 10] and dijet production in the final state 1111 
12 1 quite precise diffractive parton distributions functions (DPDF) have been extracted from 
data, which parametrize the parton content of the color singlet exchanged in the f -channel. 
The comparison of QCD predictions for single diffrative hard processes based on diffractive 
parton distributions measured at HERA and assuming factorization, against data measured 
at Tevatron [ 19.20], have indeed revealed that these processes are, not unexpectedly 1171 
[T8l . significantly suppressed in hadronic collisions, see the very recent analysis reported in 
Ref. [13]. Recalling that these distributions are fracture functions in the z — > 1 limit, the 
present formalism can then applied to next-to-leading order accuracy, the main contribution 
to the cross-sections coming from the target term, da H T . Such term can be eventually recast 
in triple differential form in xjp ~ 1 — z, virtual photon rapidity, y, and invariant mass Q 2 and 
evaluated at next-to-leading order by using the appropriate coefficient functions 11141 . In this 
way factorization tests could be porformed at fixed xip to avoid any Regge factorization 
assumption on DPDF while the y dependence, giving direct access to the fractional parton 
momentum in the diffractive exchange, jS, allows to test factorization in a kinematic region 
which avoids DPDF extrapolation. Finally, the Q 2 dependence of the cross-sections could 
be used to investigate how factorization breaking effects eventually evolve with the hardness 
of the probe and to which extent the factorized formula M ® f actually works, as anticipated 
in (2fl . 

The calculation has been performed to make predictions for cross-sections integrated 
over partonic transverse momentum. To this end, divergent contributions, related to parton 
emissions at vaninshing transverse momentum, are factorized into fracture functions. In the 
case, however, that cross-sections are measured down to a minimum but still perturbative 
hadronic transverse momentum, the latter constitutes a natural infrared regulator for the 
partonic matrix elements. The central production term, da H,c ^ l \ can be used to estimate 
hadron production as the fragmentation process, parametrized by fragmentation functions, 
were happening in the QCD vacuum. The hadronic cross-sections can be recast in a triple 
differential form in Q 2 , produced hadron transverse momentum p, and pseudo-rapidity r\ 
to predict charged particle spectra or multiplicity. A particular interesting observable which 
can also be reconstructed is the differential cross-sections differential in cos0, where (j> 
is the angle formed by the virtual photon and the detected hadron in the center of mass 
system. This observable has been shown to be sensitive to the contamination of the so-called 
underlying event 1 15 ] to jet observable and has been used also to investigate underlying event 
properties in Drell-Yan process 1161 . Although predictions made on the present formalism 
take into account the radiation accompanying one single hard scattering per proton-proton 
interactions, it can be nevertheless used as a reference cross-sections to gauge the impact of 
new phenomena, like double parton scattering. The da H C term, altough formally ff{a s ), is a 
tree level predictions and possible large higher order corrections may be expected especially 
in the forward region at large transverse momentum. In view of this fact a &{(X 2 ) calculation 
at finite transverse momentum, as performed in DIS 1221 . would be highly desiderable. 
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4 Conclusions 

We have briefly reviewed a perturbative approach to single particle production associated 
with a Drell-Yan pair in hadronic collisions. On the theoretical side we have shown that the 
introduction of new non-perturbative distributions, fracture functions, allows a consistent 
factorization of new class of collinear singularities stemming for configurations in which the 
parent parton of the observed hadron is collinear to the incoming parton. The scale depen- 
dence induced by this generalized factorization is driven by Altarelli-Parisi inhomogeneous 
evolution equations for fracture functions which allow the resummation to all orders of this 
new class of collinear logarithms. The factorization procedure does coincide with the one 
used in DIS confirming, as expected, the universal structure of collinear singularities among 
different hadron initiated processes and supporting the collinear factorization formula pro- 
posed in eq. (4). On the phenomenological side we have briefly discussed a few applications 
in which different aspects of the formalism could be tested. The improved theoretical control 
on the perturbative component indeed allows the investigation on new phenomena appearing 
in hadronic collisions, for example the rapidity gap probability suppression in hard diffrac- 
tive processes respect to diffractive DIS and the investigation, although indirect, of multiple 
parton-parton interactions. 
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